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ABSTRACT 


This  report  is  one  of  a  series  on  the  geometrical  and 
topological  properties  of  three-dimensional  digital  images. 

In  this  report,  we  define  the  genus  of  a  subset  of  a  three- 
dimensional  array  and  show  that  the  properties  of  this 
digital  genus  are  consistent  with  those  of  the  genus  in 
ordinary  (continuous)  topology.  Since  the  approach  here  is 
different  than  that  normally  used  in  two  dimensions,  we  also 
present  the  two-dimensional  case.  The  derivation  of  the 
three-dimensional  genus  allows  us  to  make  explicit  what  is 
meant  by  a  hole  in  an  object  (e.g.,  the  hole  in  a  ring) . 

We  also  give  algorithms  for  computing  the  three-dimensional 
digital  genus  and  discuss  the  importance  of  three-dimensional 
genus  in  contrast  with  the  two-dimensional  genus. 
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1.  Introduction 


The  analysis  of  three-dimensional  (3D)  imaaes  has  become 
of  increasing  interest  (see,  for  example,  (1-6])  with  the 
rapid  growth  of  computed  tomography,  in  which  discrete  3D 
representations  of  solid  objects  are  reconstructed  from  sets 
of  projections.  This  report,  in  which  we  define  (and  give 
algorithms  for  computing)  the  genus  of  a  binary  3D  image 
( as  might  be  obtained  by  applying  a  threshold  to  an  image 
produced  by  computed  tomography) ,  is  one  of  a  series  (7-10) 
on  the  geometry  and  topology  of  3D  arrays. 

Geometrical  and  topological  properties  of  digital  arrays 
play  an  important  role  in  image  analysis  and  recognition. 

While  for  two-dimensional  (2D)  images  the  genus  of  a  digital 
image  is  not  a  particularly  useful  property  per  se,  it  is  of 
theoretical  interest.  For  example,  there  is  a  well  developed 
theory  [11,12]  of  the  topological  properties  (e.g.,  connected¬ 
ness,  objects,  holes)  of  2D  digital  images  which  combine  with 
various  geometric  concepts  (e.g.,  arcs,  curves)  to  form  the 
basis  of  a  consistent  representation  scheme  that  underlies 
many  image  analysis  tasks  such  as  object  extraction  and  detec¬ 
tion,  thinning,  and  skeletonization. 

We  assume  some  familiarity  with  the  notions  of  connected¬ 
ness,  objects,  holes,  etc.,  as  in  [11,12]  for  2D  digital 
images.  Connectedness  is  a  relation  defined  on  binary  valued 
points  of  a  digital  array  which  partitions  the  array  into 


objects  (connected  components  of  points  with  image  value 
1)  and  non-objects  (connected  components  of  points  with 
image  value  0),  of  which  all  but  a  single  (unique)  background 
component  are  holes.  Simple  algorithms  may  be  given  for 
labeling  (and  counting)  objects  and  holes. 

The  genus  of  a  2D  digital  image  may  be  defined  to  be  the 
number  of  objects  minus  the  number  of  holes  in  the  image. 

While  counting  objects  only  or  holes  only  is  necessarily  a 
global  operation  [13]  there  exist  counting  functions  of  local 
patterns  of  points  which  determine  the  genus  of  the  image  [11- 
15] .  When  4 -connectedness  (see  below.  Section  3)  is  used  for 
the  l's  (a  subset  S)  of  the  image,  we  have 

G4(S)  =  *^2+*  3  (1) 

where 

ipj.  =  number  of  l's  in  the  image 

4>2  =  number  of  1x2  blocks  of  l's  in  the  image  (in  all  orien¬ 

tations) 

-  number  of  2x2  blocks  of  l's  in  the  image 
and  when  8-connectedness  is  used  for  S  we  have 

G8(S)  "  c*’l"(cf,2+93)+?4"v5  (2) 

where 

cp .  =  number  of  l's  in  the  image 
1  1 

=  number  of  times  the  pattern  11  or  1  occurs  in  the  image 

9 2  *  number  of  times  the  pattern  1^  or^l  occurs  in  the  image 

.  1111  11 

9 4  *  number  of  times  the  pattern  11,  x,  1  1,  or  1 

occurs  in  the  image 


1  1 

<Pg  *  number  of  times  the  pattern  1  1  occurs  in  the  image 

For  3D  digital  images  we  may  define  connectedness  (see 
below,  Section  4)  as  a  relation  on  binary  valued  points  of  a 
digital  array  in  a  fashion  analogous  to  the  2D  definition, 
except  that  we  call  all  but  the  background  components  of  0's 
cavities  rather  than  holes  [7,8].  The  reason  for  this  change 
in  nomenclature  will  soon  become  evident. 

In  [15]  it  was  argued  that  for  6-connected  sets  of  l's 
(see  below)  in  a  3D  digital  image  the  analog  of  Eq.  1  is 
Gg (S)  * 

where 

4^  *  number  of  l's  in  the  image 

i|»2  *  number  of  lxlx2  blocks  of  l's  in  the  image  (in  all  orien¬ 
tations) 

ij>3  *  number  of  1X2X2  blocks  of  l's  in  the  image  (in  all  orien¬ 
tations) 

^  =  number  of  2X2X2  blocks  of  l's  in  the  image 
and  (as  we  shall  see)  the  analog  of  Eg.  2  for  26-connected  sets  is. 

G26  (S>  =  *rV*3"V*5"V*7_*8 

where  the  are  as  defined  in  Appendix  A. 

Consider  the  3D  array  E^  (x,y,z)  constructed  from  a  2D 
(2) 

array  Z'  *  (x,y)  by 

E  ^  (x,y)  if  z=zQ 


I(3)  (x,y,z) 


0 


otherwise 


(2) 

where  I  has  O  objects  and  H  holes.  It  is  easy  to  see  that 

(21 

with  and  Gg  defined  on  E '  ,  and  Gg  and  G^g  defined  on 

E^,  we  have  G4=Gg  and  Gg=G2g'  and  that  the  image  E^  has  no 
cavities.  Thus,  although  cavities  are  the  3D  analog  of  2D 
holes  as  defined  by  the  connectedness  relations,  the  functions 
Gg  and  G2g  do  not  give  the  number  of  3D  objects  minus  the  num¬ 
ber  of  3D  cavities.  In  the  simple  construction  above,  we 
see  that  Gg  and  G26  are  also  affected  by  what  are  intuitively 
3D  holes,  i.e.,  holes  in  a  sheet. 

In  3D  the  genus  is  sometimes  defined  [14,15]  to  be  the 
number  of  objects  minus  the  number  of  holes  plus  the  number  of 
cavities.  Since  for  3D  digital  images  the  topological  con¬ 
cepts  of  object  and  cavity  are  well  understood  (in  particular, 
algorithms  similar  to  those  for  2D  images  can  be  given  for  count¬ 
ing  them)  a  definition  of  genus  would  define  the  number  of 
holes,  and  conversely.  What  we  desire  is  some  evidence  to 
support  the  assertion  of  either  definition  as  being  consistent 
with  the  corresponding  ordinary  topological  notion. 

Even  in  ordinary  topology  it  is  difficult  to  characterize 
holes.  A  hole  may  be  thought  of  as  a  property  of  a  boundary 
surface  (separating  binary  regions)  which  makes  it  topologically 
equivalent  to  a  torus;  an  object  with  H  holes  has  a  surface 
topologically  equivalent  to  an  H-holes  torus.  In  another 
approach,  an  object  is  defined  to  have  no  holes  if  every 
simple  closed  curve  in  the  object  is  continuously  deformable 
within  the  object  to  a  single  point  (see  [7]  where  this  approach 


is  used  for  digital  images;  for  the  present  purposes  we 
need  a  stronger  statement/  since  we  will  in  general  be  deal¬ 
ing  with  objects  that  may  have  many  holes) . 

We  see  from  these  remarks  that  the  concept  of  a  hole  is 
different  from  those  of  objects  and  cavities;  we  cannot 
point  to  or  label  the  points  which  constitute  a  hole.  Indeed, 
the  points  of  objects  and  cavities  cover  the  space,  so  that 
a  hole  is  a  property  of  these  collections  of  points.  Thus, 
when  considering  an  object  (and  its  cavities)  we  shall  here 
try  only  to  understand  what  is  meant  by  the  number  of  holes 
in  the  object,  and  not  what  is  meant  by  a  hole. 

We  begin  by  postulating  the  well-known  Euler-Schlaef li 
relation  as  a  fundamentally  true  statement  of  dimensional 
connectedness.  This  equation  describes  a  relationship  between 
the  numbers  of  elements  of  different  dimensionality  in  an 
interlinked  system  of  these  elements  given  the  geometrical 
constraints  of  Euclidean  space.  This  theory  of  interlinked 
systems  of  elements  (structures)  requires  that  these  structures 
have,  in  effect,  no  holes  (i.e.,  are  singly  connected).  We 
have  chosen  this  starting  point  because  the  process  of  extend¬ 
ing  the  concepts  of  singly  connected  structures  to  more  complex 
(multiply  connected)  structures  furthers  an  understanding  of 
what  is  meant  by  the  definition  of  digital  genus  (number  of 
holes) ,  and  not  because  multiply  connected  structures  are  not 


also  well  understood. 


The  approach  taken  here  is  constructive.  That  is,  a 
structure  is  described  by  specifying  the  elements (of  differ¬ 
ent  dimensionalities)  that  constitute  it.  We  proceed  by 
making  two  associations:  one  from  points  of  an  image  to 
structures,  and  one  from  components  (objects,  holes (2D) ,  cavi¬ 
ties  (3D)  )  to  elements  of  the  structures.  Even  for  2D  digital 
images,  this  approach  is  essentially  different  than  the  expo¬ 
sition  of  genus  found  in  [11,12,13].  Thus,  following  the 
introduction  of  structures  and  the  Euler-Schlaef li  identity 
(Section  2)  we  illustrate  this  technique  by  presenting  a  deri¬ 
vation  of  genus  for  2D  digital  images  (Section  3)  before 
developing  the  genus  (number  of  holes)  of  a  3D  digital  image 
(Section  4).  In  Section  5,  we  give  an  algorithm  for  computing 
the  3D  digital  genus  of  an  image. 

Some  early  work  on  topology  of  3D  digital  arrays  is 
found  in  [14,15],  and  several  theoretical  papers  [16,17]  also 
consider  generalizations  to  higher  dimensions.  More  recently 
[18-20]  discuss  the  topology  of  surfaces  of  3D  objects. 


2 .  Space  Structures 

The  theory  of  structures  given  here  is  from  the  excellent 
book  by  Loeb  [21] .  In  general  a  structure  is  an  underlying 
pattern  of  interaction  between  things.  When  a  structure  is 
confined  to  a  Euclidean  space,  geometrical  constraints  are 
imposed  on  the  structure,  which  allow  us  to  write  various 
quantitative  expressions  of  spatial  order.  The  well  known 
Euler-Schlaef li  equation 

j  i  i 

E  (-1)  N.  =  1  +  (-l)3 

i=0 

describes  the  relationship  between  the  numbers  of  elements  of 
different  dimensionality  in  a  multidimensional  singly  connected 
structure  confined  to  Euclidean  space.  Here  is  the  number 
of  elements  of  dimensionality  i  and  j  is  the  dimensionality 
of  the  structure.  [In  two  dimensions  we  write  this  as  the 
more  familiar  Euler  invariant  for  any  planar  graph, 

V-E+F  =  2 

where 

V  =  Nq  =  number  of  vertices 
E  =  =  number  of  edges 

F  -  Nj  =  number  of  faces.] 

We  note  that  in  this  identity  the  infinite  region  surrounding 
the  structure  is  considered  an  element  of  dimensionality  N ^ . 

Some  remarks  about  the  use  of  this  identity  are  necessary. 
First,  dimensionality  refers  to  numbers  of  degrees  of  freedom. 


A  point  confined  to  move  inside  a  cube  or  sphere  may  move 
in  three  mutually  perpendicular  independent  directions.  A 
point  confined  to  the  surface  of  a  cube  or  of  a  sphere  (even 
though  not  flat)  has  only  two  independent  directions  of  motion. 
Similarly,  a  line  or  curve  has  dimensionality  1  and  a  point 
has  dimensionality  zero. 

Next,  this  identity  applies  only  to  structures  that  are 
singly  connected.  A  structure  is  multiply  connected  if 
there  are  two  or  more  curves  between  any  two  points?  which 
cannot  be  brought  into  alignment  by  continuous  deformation 
(cf.  the  definition  of  the  class  of  objects  with  no  holes 
in  [7]).  This  has  a  profound  effect  on  the  way  we  must  view 
structures.  Consider,  for  example,  the  2D  structure 

A 

The  singly  connected  restriction  forces  us  to  view  the  region 
surrounded  by  the  vertices  and  edges  as  a  face  of  the  structure 
(the  exterior  is  also  considered  a  face) .  Just  as  an  edge  is 
commonly  used  to  represent  an  independent  relation  between  two 
vertices  in  a  graph,  the  face  in  the  above  structure  represents  an 
independent  relation  between  each  vertex  and  the  other  two. 

In  the  3D  structure  (tetrahedron) 


we  are  left  a  choice:  The  (3D)  structure  consists  of  4  vertices, 
6  edges  together  with  either  3  faces  and  1  "background"  cell 


(structured  element  of  dimensionality  3) ,  so  that  the  struc¬ 
ture  is  cup-like,  or  4  faces  so  as  to  create  a  surface  that 
divides  the  three-dimensional  space  into  an  enclosed  cell 
and  all  the  space  outside  it. 

Additionally,  we  shall  require  that  the  following  closure 
postulate  be  satisfied:  at  least  two  edges  join  at  every 
vertex,  and  at  least  two  faces  meet  at  every  edge.  Thus  we 
exclude  tree  graphs  from  consideration  as  structures  (although 
it  is  possible  by  adopting  a  peculiar  vantage  point  to  view 
tree  graphs  as  structures  [21,pp. 8 , 9] ) ,  and  in  structures 
like  that  above  the  choice  of  structures  is  removed.  Never¬ 
theless,  we  are  left  with  a  choice  in  structures  such  as 


Either  there  are  ten  faces  and  two  cells  or  eleven  faces  and 
three  cells. 


3.  2D  Digital  Images 

A  2D  digital  image  S  is  a  2D  lattice  of  points  defined 
by  pairs  of  Cartesian  coordinates  (x,y)  which  we  may  take  to 
be  integer  valued.  A  binary  2D  digital  image  is  a  2D  digital 
image  whose  points  have  values  0  or  1.  The  geometrical  and 
topological  concepts  presented  here  may  be  found  in  [11,12]. 

We  will  consider  two  types  of  neighbors  of  a  point  p=(x,y): 

(a)  the  neighbors  (u,v)  such  that  |x-u|  +  |y-v|  =  1, 
which  we  call  edge  neighbors  of  p; 

(b)  the  neighbors  (u,v)  such  that  |x-u|  =  1  and  |y-v|  =  1, 
which  we  call  corner  neighbors  of  p. 

The  four  edge  neighbors  of  p  are  said  to  be  4-adjacent  to  p, 
and  the  eight  edge  and  corner  neighbors  of  1  are  said  to  be 
8-ad jacent  to  p. 

A  path  from  p  to  q  is  a  sequence  of  points  p=pQ , . . . »Pn=q 
of  E  such  that  pi  is  adjacent  to  p^_1#  1  *  i  *  n,  and  any  point 
alone  is  a  path  of  length  zero.  In  general,  for  any  r on-empty 
subset  S  of  E  we  speak  of  (4  or  8) -paths  in  S  depending  on  the 
type  of  adjacency  used. 

Let  S  be  a  subset  of  points  of  Z  and  to  avoid  special  cases 
assume  that  S  does  not  meet  the  border  of  Z  .  We  say  that  p  and 
q  are  connected  in  S  if  there  is  a  path  from  p  to  q  in  S.  Readily, 
connectivity  is  an  equivalence  relation,  and  the  equivalence  clas¬ 
ses  under  this  relation  are  called  components  (or  objects)  of  S. 


Note  that  these  are  each  two  definitions  in  one,  depending  on  the 
type  of  adjacency  used,  4-  or  8-. 


Similarly  we  can  consider  components  of  the  complement 
S  of  S.  Evidently  exactly  one  of  these  components  contains  the 
border  of  E;  we  call  this  component  the  background  of  S.  All 
other  components  of  S  (if  any)  are  called  holes  in  S.  To  avoid 
ambiguous  situations,  we  will  require  that  opposite  types  of 
connectedness  be  used  for  S  and  S.  Simple  algorithms  have  been 
given  to  label  and  count  components  of  S  and  S  [11]. 

L^t  A, B  be  subsets  of  I.  We  say  that  A  surrounds  B  if 
any  path  from  (a  point  of)  B  to  (a  point  of)  the  border  of 
E  must  meet  (contain  a  point  of)  A.  From  [12]  we  have 
Theorem  1.  Any  S  surrounds  its  holes  and  is  surrounded  by  its 
background. 

Theorem  2.  The  adjacency  graph  of  S  (the  graph  whose  nodes 
are  components  of  S  and  S  and  whose  edges  are  defined  by  adjacency 
of  components  of  S  and  Sj  is  a  directed  tree  rooted  at  the  back¬ 
ground  component  of  S  under  the  relation  of  surroundedness. 

These  two  facts  are  significant  in  developing  an  algorithm 
to  compute  the  genus  of  a  set  S  (or  S)  of  a  digital  image. 

Namely,  we  can  count  the  components  of  S  as  objects  and 
count  each  set  that  is  surrounded  by  a  component  of  S  as  a 
hole  to  compute  the  genus.  Thus  it  suffices  to  consider  the 
genus  of  a  single  component  of  S  and  its  complement. 

With  every  point  of  a  digital  image  we  may  associate  a  pixel 
(or  "cell") ,  which  is  a  square  region  defined  by  four  corners 
and  four  edges  which  are  shared  with  the  pixels  of  neighboring 
points.  We  associate  the  binary  value  of  each  point  with  the 


face  of  tae  pixel. 


The  genus  G(S)  of  a  set  is  the  number  o(S)  of  objects 
(components  of  S)  minus  the  number  H(S)  of  holes  in  S.  Thus, 
the  component  labeling  and  counting  algorithms  may  be  used  to 
determine  tne  genus  of  a  set  S.  While  the  labeling/counting 
algorithms  are  global  and  sequential  in  nature,  they  have 
been  used  to  demonstrate  the  existence  of  local  parallel 
counting  functions  that  determine  the  genus  by  an  induction 
on  the  number  of  points  in  S[12];  given  any  n  point  object 
for  which  the  local  counting  functions  agree  with  the  global 
counting  algorithms,  it  can  be  shown  that  adding  a  point  to 
the  image  will  have  the  same  effect  on  both. 

We  will  see  later  that  this  approach  cannot  be  used  for 
3D  digital  images.  In  3D  we  can  present  algorithms  to  label 
and  count  objects  (components  of  S)  and  cavities  (non¬ 
background  components  of  S) ,  but  not  holes;  a  hole  is  not  some¬ 
thing  to  which  we  can  assign  a  label  in  the  same  sense  that 
we  label  objects  or  cavities. 

Consider  the  genus  G(S)  of  the  complement  S  of  S  in  a  2D 
digital  image.  We  see  (from  the  adjacency  tree)  that  every 
object  in  S  is  a  hole  in  S  (i.e.,  0(S)  =  H(S))  and  that  every 
object  of  S  except  the  background  is  a  hole  in  S  (i.e.,  0(5)-l 
=  H (S ) ) .  We  thus  have 

0  (S )  -  H  (S)  =  H(S)  -  0(S)  +  1 
or 

G  (S)  +  G(S)  =  1 

Since  we  use  opposite  types  of  connectedness  for  S  and  S,  we 
see  that  the  digital  genus  we  define  for  each  type  of  connected¬ 
ness  must  satisfy  G.  (S)  +  G_(S)  =  1  and  GQ(S)  +  G.  (S)  = 

4  o  o  4 


i. 


We  first  consider  the  genus  of  a  single  8-component 
C&S.  Consider  the  graph  whose  vertices  are  corners  of  pixels 
in  S  and  whose  edges  are  edges  of  these  pixels,  e.g.,  from 
the  image 

111 
1  11 
111 

we  obtain  the  graph 


lot 

l 

a 

a" 

■a  16 
( 

a 

a 

a 
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(Clearly  every  such  graph  is  connected  and  planar.)  Viewing 
the  graph  as  a  structure  we  can  identify  (specify)  three 
distinct  types  of  faces: 

a)  faces  arising  from  pixels  in  C 

3)  faces  arising  from  components  of  £  surrounded  by  C 
)  the  "background"  face  arising  from  the  component  of  C 
not  surrounded  by  C. 

Clearly,  every  pixel  of  C  becomes  a  face  of  the  graph  by 
construction.  The  single  component  of  C  which  surrounds  C  is 
associated  with  the  "background"  face  of  the  structure.  Every 
other  component  Due  is  surrounded  by  C.  Define  an  adjacency 
relation  between  the  component  D  and  the  pixels  of  the  D-border 
of  C,  and  we  see  that  we  may  associate  each  component  D  with  the  face 
of  the  structure  adjacent  to  the  faces  arising  from  the  pixels  of 
the  D-border  of  C.  The  face  associated  with  each  D  is  unique  since 
the  D-border  of  C  surrounds  D.  In  general  faces  of  type  (6)  are 


n-gons  with  n<s  4.  We  can  new  write 
V-E+F  -  V-E+(#a  +  #fl  ♦  #Y ' 


2 


or 

1-#B  =  V-E+#a 

where  #a,  #S,  #y  are  the  numbers  of  these  types  of  faces 
in  the  structure. 

Now  #B  is  the  number  of  faces  resulting  from  components 
of  C  surrounded  by  C,  i.e.,  holes  in  C.  As  we  have  seen, 
it  is  sufficient  to  sum  the  number  of  holes  in  each  C&S  to 
determine  the  number  of  holes  in  S.  Thus,  if  we  sum  V,  E, 
and  #a  over  the  structure  of  every  component  C&S  we  have 

0 (S) -H (S)  =  V-E+#a 

where  0(S)  is  the  number  of  objects  and  H(S)  is  the  number  of 
holes  in  S. 

We  now  consider  computing  V-E+#a  directly  from  the  digi¬ 
tal  image  I.  Clearly  #a  is  just  the  number  of  pixels  with 
image  value  1.  Similarly,  it  is  not  hard  to  see  that 

V  »  4<f.  -  2Y0  -  +  <fr,  -  Cfc 

1  2  3  4  5 

E  «  -  *2 

so  that 

V-E+#a  =  -  ('*2+'r3)  +  V4  -  v 5 

where 

y,  =  number  of  points  in  S 

1 

*2  =  number  of  times  the  pattern  1  1  or  1  occurs  in  the  image 

y3  =  number  of  times  the  pattern  l^-  or  ^1  occurs  in  the  image 

,  .  1  1  1  111 

*  number  of  times  the  pattern  1  1  or  1  or  1  1  or  1 

occurs  in  the  image 


1  1 

=  number  of  times  the  pattern  1  1  occurs  in  the  image. 

We  thus  have 

Gg(S)  =  0(S)  -  H(S)  =  <P1-('P2+'P3)-»-94-'P5 

as  the  genus  of  S  when  we  use  8 -connectivity  for  S. 

Several  approaches  suggest  themselves  for  dealing  with  the 
4-connected  case.  For  example,  we  have  already  noted  that 
G(S)+G(S)  =  1.  Consider  the  structure  of  E  as  per  the  construc¬ 
tion  above.  It  is  easy  to  see  that  for  an  mxn  array  of  points 
we  count  V' ,  E' ,  and  F*  of  this  structure  as 
V'  =  (n+1) (m+1) 

E '  *  2mn+m+n 
F '  =  nm+1 
so  that 

V'-E'+F'  «  2 

Counting  directly  from  the  image  as  before  we  can  write 
V'  =  4<P1-2cfr2-93+5P4-V5+v3+2m+2n 
E'  =  4<f>1-V2+iJ>2+2m+2n 
F'  * 
so  that 

V’-E'+F*  =  ['f'L-(^2+S'3)+cf,4-^5]+[i|/1-tp2+^3]+l  =  2 
where 

ip  =  nximber  of  points  in  S 

0 

4>2  -  number  of  times  the  pattern  0  or  0  0  occurs  in  the  image 

Q  0 

t^3  =  number  of  times  the  pattern  o  0  occurs  in  the  image 
Thus,  we  have  Gg  (S)  +  [^1-^2+i};3)  =  1,  so  that  G4  (S)  *^^-^2+^2 .  We 
shall  see  that  for  3D  digital  images  this  approach  requires 
an  additional  assumption. 


Thus,  let  us  proceed  as  before  by  constructing  a  structure 

which  corresponds  to  a  4 -component  C  of  S.  The  construction 

0  1 

is  similar  except  that  for  the  configuration  1  0  (and  its 
reflection)  we  construct  the  subgraph 


i.e.,  there  are  two  vertices  at  the  center  of  the  configuration. 

The  proof  that  we  can  identify  three  distinct  types  of  faces  of 
the  structure  (ot-y)  stands  without  change,  and  we  obtain 
O(S)  -  H (S)  =  V-E+#a 

We  can  count  V  and  E  directly  from  the  image  as 
v  =  4iJ*1  -  24>2  +  ip3 
E  =  44>l  -  ip2 
#a  ■ 
so  that 

V-E+#a  =  4>1-ip2+^3 

where  the  are  as  before  but  for  image  values  of  1.  We  thus  have 
G4(S)  =  0  (S)  -H  (S)  =  ^-^2+4^ 

Even  with  the  complication  of  the  additional  vertex  for  a  pair 
of  diagonal  l's  the  calculations  of  the  preceding  paragraph  stand 
without  change,  so  that  we  could  in  principle  define  the  genus 
Gg (S)  of  the  complement  in  terms  of  G^ (S)  by  considering  the 
structure  corresponding  to  E. 

Before  discussing  the  genus  of  a  3D  digital  image  in  the  next 


section  we  should  point  out  that  the  structures  used  here  for  2D 
digital  images  are  by  no  means  the  only  structures  we  can  use 


for  defining  genus.  For  example,  it  is  not  hard  to  see  that 


*  4^-41 2+ij>3  may  be  interpreted  as  counting  vertices,  edges, 
and  faces  on  a  graph  whose  vertices  are  points  of  the  image 
and  whose  edges  are  defined  by  4-connectivity,  although  a  some¬ 
what  different  association  must  be  made  between  faces  of  the 
graph  and  subsets  of  the  digital  image.  Note  that  when  the  edges 
of  such  a  graph  are  defined  by  8-connectivity  the  "structure" 
becomes  three-dimensional.  For  3D  digital  images  this  kind  of 
graph  will  lead  to  an  efficient  algorithm  for  computing  the  genus. 


L 


1 


.  Digital  Binary  Images  in  Three  Dimensions 


A  three-dimensional  digital  image  Z  is  a  three-dimensional 
lattice  of  points  defined  by  triples  of  Cartesian  coordinates 
(x,y,z)  which  we  may  take  to  be  integer  valued.  A  binary  3D 
digital  image  is  a  digital  image  whose  points  have  values  0  or 
1.  In  this  section  we  review  the  concepts  of  adjacency,  connected¬ 
ness,  and  components  of  points  in  Z  as  developed  in  [7,8,10]. 

We  will  consider  two  types  of  neighbors  of  a  point  p  = 

(x,y,z) : 

(a)  the  neighbors  (u,v,w)  such  that  j x-u | + | y-v | + | z-w |  =  1, 
which  we  call  6-neighbors  of  p 

(b)  the  neighbors  (u,v,w)  such  that  max[|x-u|,  | y— v | ,  | z-w | ]  =  1, 

which  we  call  the  2 6 -neighbors  of  p. 

The  6-neighbors  are  said  to  be  6-adjacent  to  p,  and  the  26-neighbors 
26-adjacent  to  p. 

A  path  from  p  to  q  in  Z  is  a  sequence  of  points  p=pQ , .  .  . ,  Pn=<3 
of  Z  such  that  p^  is  adjacent  to  p^_^,  l*i*n.  Any  point  alone 
is  a  path  of  length  zero.  In  general,  for  any  non-empty  subset 
S  of  Z  we  speak  of  (6-  or  26-)  paths  in  S  depending  on  the  type 
of  adjacency  used. 

Let  S  be  a  subset  of  points  of  Z ,  and  to  avoid  special  cases 
assume  that  S  does  not  meet  the  border  of  Z .  We  say  that  p  and 
q  are  connected  in  S  if  there  is  a  path  from  p  to  q  in  S.  Readily, 
connectivity  is  an  equivalence  relation,  and  the  equivalence  clas¬ 
ses  under  this  relation  are  called  components  of  S  (or  objects) . 

Note  that  these  are  two  definitions  in  one,  depending  on  the  type 
of  adjacency  used. 


Similarly;  we  can  consider  components  of  the  complement 
S  of  S.  Evidently  exactly  one  of  these  components  contains  the 
border  of  E;  we  call  this  component  the  background  of  S.  All 
other  components  of  S  (if  any)  are  called  cavities  in  S.  To 
avoid  ambiguous  situations  we  will  require  that  opposite  types 
of  adjacency  be  used  for  S  and  S.  Simple  algorithms  can  be  given 
for  labeling  and  counting  distinct  components  of  S  and  S  [15,7]. 

The  genus  G(S)  of  a  set  S  in  a  3D  digital  image  is  the  number 
0 (S)  of  objects  plus  the  number  C(S)  of  cavities  in  S  minus  the 
number  H(S)  of  holes  in  S.  As  already  mentioned  the  definition 
of  hole  is  not  simple,  and  in  particular  holes  cannot  be  labeled 
to  facilitate  counting  them,  as  can  be  done  with  objects  and 
cavities. 

Let  A,  B  be  subsets  of  E.  We  say  that  A  surrounds  B  if  any 
path  from  (a  point  of)  B  to  (a  point  of)  the  border  of  E  must  meet 
(contain  a  point  of)  A.  As  in  2D  we  have  [  7  ] 

Theorem  3.  Any  S  surrounds  its  cavities  and  is  surrounded  by 
its  background. 

Theorem  4.  The  adjacency  graph  of  S  (the  graph  whose  nodes  are 
components  of  S  and  S  and  whose  edges  are  defined  by  adjacency 
of  components  of  S  and  S)  is  a  directed  tree  rooted  at  the  back¬ 
ground  component  of  S  under  the  relation  of  surroundedness. 

With  every  pont  of  a  digital  image  we  may  associate  a  voxel . 
usually  a  cubic  region  defined  by  8  corners,  twelve  edges,  and 
six  faces  which  are  shared  with  the  voxels  of  neighboring  points. 
We  thus  associate  the  binary  value  of  each  point  with  the  cell 
(or  volume)  of  the  voxel. 


For  any  3D  structure  the  Euler-Schlaef li  equation  is 
V-E+F-Q  =  0 
where 

V  =  number  of  vertices 

E  =  number  of  edges 

F  =  number  of  faces 

Q  =  number  of  cells  (including  the  background) . 

From  a  3D  digital  image  we  may  construct  a  structure  from  any 
single  27-component  Bcs  whose  vertices  are  corners  of  voxels 
with  image  value  1,  whose  edges  are  edges  of  these  voxels,  and 
whose  faces  are  faces  of  these  voxels.  In  analogy  with  the  2D 
development  we  can  identify  three  types  of  cells: 

(a')  cells  arising  from  voxels  of  B 

(3')  cells  arising  from  components  of  B  surrounded  by  B 
(<5')  the  "background"  cell  arising  from  the  component  of  B 
not  surrounded  by  B. 

In  general,  cells  of  type  3'  are  polyhedrons  with  eight  or  more 
vertices.  The  details  of  this  association  are  similar  to  those 
in  the  2D  proof. 

This  done,  we  find  that  the  resulting  structure  is  not  neces¬ 
sarily  singly  connected.  As  a  simple  example  consider  the  structure 
resulting  from  the  digital  image 

0  110 
10  0  1 
0  110 


which  is 


Some  additional  structural  elements  are  required  to  make  this 
structure  singly  connected  (cf.  the  "missing  faces"  of  the 
2D  development  above) . 

In  general,  for  the  structure  derived  as  above  from  any 
digital  component  BtS  we  may  compute  V-E+F  (#a+#3+#y)  although 
this  sum  need  not  equal  zero  since  the  structure  need  not  be 
singly  connected.  We  define  the  number  of  holes  in  the  object 
as  the  amount  by  which  this  sum  differs  from  zero, 

H (B)  =  -V  +  E  -  F  +(#a+#S+#Y) 

Noting  that  #0  is  the  number  of  cavities  C(B)  in  any  component 
BcS  and  that  #y  =  1  (for  the  background)  we  have 
1+C (B) -H (B)  =  V-E+F-#a 

so  that  we  may  sum  V,E,F,#a  over  every  component  B&S  in  the  image 
to  obtain 

0 (S) -H (S) +C (S)  =  V-E+F-#a 

where  0(S)  is  the  number  of  objects  in  the  image,  C(S)  is  the 
total  number  of  cavities  in  these  objects,  and  H(S)  is  the  total 
number  of  holes  in  these  objects.  [We  are  tacitly  assuming  that 
H (S)  =  EH (B) ,  where  the  summation  is  over  the  components  B  of  S. 

This  additive  property  of  holes  of  components  is  justified  by  the 
fact  [14]  that  genus  is  a  local  set  property,  viz.  G(A)+G(B)  = 
G(AHB)+G(AUB)  for  objects  A  and  B.] 

Remark:  The  definition  of  number  of  holes  above  may  be  related 

to  the  Euler-Poincare  characteristic  for  finite  complexes  [22] . 

The  structure  corresponding  to  a  digital  component  is  a  complex 
K-L  where  L  is  a  subcomplex  of  a  simplicial  complex  K.  Intuitively, 


the  sum  V-E+F- ( #a+#8+#y)  differs  from  zero  by  some  amount 
v-e+f-g,  where  these  are  computed  on  K-L  and  L,  respectively. 
Viewed  in  this  fashion  the  number  of  holes  is  the  sum  v-e+f-g 
for  any  structural  elements  added  to  the  structure  resulting 
from  a  digital  image  which  will  make  it  singly  connected. 

We  now  consider  computing  V-E+F-#a  directly  from  the  digi¬ 
tal  image.  Although  more  laborious  than  the  2D  case,  we  can 
show  that 


G27  (S)  =0  (S)  -H  (S)  +C  (S)  =c?1-Cf2  +  CVCf4+Cp5-CP6+CP7“Cp8 


where 


4'1  =  #[1] 

*2  =  #  [2]+#  [3]+# [4] 
=  #  [5]+#[6]+#[7] 


=  #l8]+#[9]+#[10]+#[ll]+#[12]+#[13]+#[14] 


*5  =  # [15]+# [16]+#  [17] 

96  =  # [18]+# [19]+#  [20] 

=  #[21] 

=  #[22] 

and  by  # [n]  we  mean  the  number  of  times  the  configuration  n  of 
Appendix  A  occurs  in  the  picture  (in  all  orientations) . 

As  in  the  2D  case  one  way  of  proceeding  to  the  6-connected 
case  is  to  note  that  for  the  structure  corresponding  to  Z  we  have 
V'-E'+F'-Q'  =  0,  and  that  counting  directly  from  the  image  we  have 
V-e’+F’-Q'  =  (91-cP2+93-cP4  +  95-96+cP7-cPg)-^1+^2-^3+^4  +  l  =  0 


(^-^2+^-1^) -G27(S)  =  1 


or 


where  the  count  patterns  of  0's  as  follows: 

=  #111 

i>2  =  #  [2] 

^3  =  #18] 

=  #[20] 

Now,  every  object  of  S  is  a  cavity  in  S  (i.e.,  0(S)=C(Sj) 
and  every  object  of  S  except  the  background  is  a  cavity  of  S 
(i.e.,  0 (S ) -1=C (5 ) ) .  If  we  make  the  assumption  that  H(S)=H(S)  we 
have  0(S)+C(S)-H(S)  =  C (S) +0 (S) -1-H (S)  or  G (S) -G (S) =1.  Thus  from 
the  above  equation  we  can  define  Gg  (S)  .  Because  we 

need  this  additional  assumption,  we  consider  it  safer  to  base  the 
definition  of  genus  for  6-connectivity  on  a  structure  correspond¬ 
ing  to  a  6-connected  object. 

The  construction  of  the  structure  for  a  6-connected  object 

is  similar  to  that  above,  except  that  for  configurations  such 
0  1 

as  1  0  and 

1st  plane  2nd  plane 

0  1  0  0 

0  0  10 

we  use  the  structures 


so  that  each  6-component  forms  a  distinct  structure.  (We  are 
glossing  over  the  details  of  this  construction;  the  essential 
point  is  that  a  vertex  is  not  shared  unless  the  voxels  are  con¬ 
nected  inside  a  2x2x2  neighborhood  containing  the  two  voxels.) 


Again,  we  identify  some  of  the  cells  of  the  structure  with 
components  of  the  complement,  and  define  the  number  of  holes  as 
the  amount  by  which  the  sum  V-E+F- (#a+#e+#y)  differs  from  zero 
(or  as  the  net  contribution  to  the  sum  V-E+F-Q  needed  to  make 
the  structure  simply  connected  i.e.,  V-E+F- (#a+#£+#Y) +H  =  0). 

For  any  single  component  BcS  we  have 
1+C (B)  -  H (B)  =  V-E+F- # a 

so  that  for  O(S)  components  of  a  digital  image  we  can  determine 
0(S)-H(S)+C(S)  =  V-E+F-#a 

by  summing  V,E,F  and  #a  over  the  entire  image.  Counting  directly 
from  the  image  we  find 
V  =  8t1-4i^2+2ilJ3-ip4 
E  =  12i^1-4^2+iJ/3 
F  =  6^1-i|»2 
#a  = 

so  that 

Gg  (S)  =0  (S)  -H  (S)  +C  (S)  =^1-^2+'^3-^4 

where  the  are  defined  as  before  but  for  patterns  of  l's  in  Z. 
Once  again,  we  find  that  (with  the  additional  assumption  H(S)=H(S)) 
we  can  in  principle  define  the  genus  of  the  complement  in  terms 
of  Gg (S)  by  considering  the  structure  corresponding  to  Z. 

From  the  definitions  of  G27  amd  Gg  based  on  structure  we  see 
that  G27(S)-Gg(S)  =  1  and  Gg(S)-G27(S)  =  1  from  counting  elements 
in  the  structure  corresponding  to  Z.  Thus,  we  have  shown  that 
H(S)=H(S)  since  0(S)=C(S)  and  0  (S)  -1=C  (S)  . 


The  equations 

Gg  (S)  =  ^1-^2+^3-if)4 

G27  (S)  ,cp1-cp2+cp3-cp4+cp5-cp6+cp7-cp8 

define  the  genus  of  a  3D  digital  image  but  are  not  algorithms 
for  computing  the  genus.  An  algorithm  for  computing  Gg(S) 
is  not  difficult  to  devise,  since  the  patterns  of  points  of  S 
Eor  ip.  and  ip4  occur  in  an  orientation  independent  fashion  and 
the  patterns  for  ^2  and  occur  in  three  different  orienta¬ 
tions  each.  Thus  it  is  sufficient  to  test  for  occurrences  of 
nine  distinct  patterns  in  the  digital  image.  For  G27(S),  how¬ 
ever,  there  are  22  distinct  patterns,  many  of  which  occur  in 
many  orientations. 

In  this  section  we  give  an  algorithm  for  computing  G27(S). 

The  development  of  this  algorithm  is  understood  in  terms  of 
a  graph  defined  by  the  27-connectivity  relation  on  the  points 
of  S  in  L  We  point  out  that  the  graph  used  here  is  not  simi¬ 
lar  to  the  structures  used  to  define  the  genus,  but  is  analo¬ 
gous  to  the  alternative  structure  for  2D  images  mentioned  at 
the  end  of  Section  3 .  The  o ^  may  be  given  the  following 
interpretation.  Consider  a  graph  <S,E(S)>  whose  vertices 
are  points  of  S  and  whose  edges  are  defined  by  the  26-connecti¬ 
vity  relation.  Then  cp^  is  the  number  of  vertices  and  cp2  the  num¬ 
ber  of  edges  of  <S ,E (S) >.  In  general,  since  every  point  in  a 

2x2x2  cube  is  26-adjacent  to  every  other  point  in  the  2x2x2 


n2=#[2] ({(i-l,j-l,k) , (i, j-lfk) , (i-l,j-l,k-l) ,<i, j-l,k-l) >) 
=  #12]  (K?  .  ) 

n3=#[2] ( { (i-1 , j , k-1) , (i,j,k-l) , (i-l,j-l,k-l) , (i, j-l,k-l) )) 


#12]  (Kf  .  .  ) 

If 


12 


n12=#[2] ({ (i-lf j-l,k) , (i-1, j-l,k-l) })  =  #[2](K*‘.  k) 
n13=#[2)  ({  (i-1, j, k-1)  ,  (i-l,j-l,k-l)})  =  #(2](kJ3.^) 

n23=#  t2]  ({  (i,  j-l,k-l)  ,  (i-1 , j-1 ,k-l)  })  =  #[2](k“.^) 

i  o  *5 

n123=#[2]  ({  (i-i/D-i.k-l) })  5 

where  # [n] (X)  refers  to  the  number  of  times  the  pattern  n  of 
Appendix  A  (in  all  orientations)  occurs  in  X. 

We  can  now  solve  for  the  net  contribution  Am„ (K.  .  ,)  of 

K.  .  ,  to  each  <*>  (S)  .  Clearly  the  1  point  objects  in  1^.  .  , 

1,J,K  X  1  ,  J  ,  K 

2 

are  counted  in  K.  ,  .  ,  the  l  point  objects  in  K.  .  .  are 
counted  in  K.  .  ,  .  and  the  £  point  objects  in  K?  .  ,  are 

l  r  3  #  k 

counted  in  K.  .  v  , .  Thus  the  count  (^)  must  be  adjusted  to 

(")-(S1)_(?2)"(?3) 


.  .  12  13 

But  now  the  numbers  of  £  point  objects  m  K.  .  .  K.  .  .  and 

i/] »»» 

23 

K.  .  ,  have  been  subtracted  twice,  so  this  must  be  further 

J 

adjusted  to 

(?)'(£1)“(?2H"3)  +  (?12)  +  (?13)  +  (S23) 

.  123 

Finally,  the  £  point  objects  m  K.  .  v  have  now  been  counted 
one  too  many  times,  so  the  final  adjustment  is 

Acp£(Ki,j,k)  =  (£)_(?1)"(?2)"(£3)  +  (£12)  +  (A13)  +  (£23)"(?123) 


This  net  contribution  of  each  2*2x2  cube  may  be  computed  for 
each  cube  of  the  image  in  parallel  and  summed. 


Similarly  this  gives  a  simple  expression  for  the  net  con¬ 
tribution  AG_C(K.  .  .  )  of  K.  .  .  to  G_C(S)  which  may  be  com- 

ZO  i  f  ]]  f  K  ±  t  J  f  K  c.  D 

puted  for  each  K.  .  .  in  parallel  and  summed: 

1  t  J  9  K 

AG26(Ki,j,k»*5  -  ‘-1) 1  <S>  -  ‘S1-)  -  ~  <?3»  + 1?12»  +  <;i3)  -M223I  - 

<"123>1 

Since  the  expression  inside  the  brackets  requires  a  significant 
amount  of  computation,  we  suggest  the  use  of  a  lookup  table 
as  follows.  From  a  2*2*2  cube  of  binary  points  form  an  8-bit 
number  t.  Clearly  the  values  of  n,  n^,  n2,  n3»  n12,  n^3,  n23, 
and  n123  may  be  predetermined  for  every  0st*255  and  the  values 
of  AG26 (<t>)  stored  in  a  table,  where  <t>  is  the  2*2*2  cube 
of  points  whose  8-bit  number  is  t.  In  Appendix  B  we  list  the 
values  of  AG2fi (<t>)  for  one  method  of  forming  the  8-bit  number 
t  from  the  8!  different  ways  of  doing  this.  It  should  be 
noted  that  in  this  table  -2*AG2g (<t>) £  1,  and  that  most  of  the 
entries  are  zeros,  so  that  if  necessary  the  size  of  the  table 
may  be  significantly  reduced. 

Of  course,  there  is  another  alternative  to  computing  the 
genus  of  a  27-connected  subset  S:  Namely,  use  the  simpler 
6-connected  algorithm  to  determine  the  genus  of  S,  and  employ 
the  identity  G(S)-G<S)=  1. 


6 .  Conclusion  and  Remarks 

We  have  proposed  definitions  of  digital  image  subsets  based 
on  the  association  of  a  structure  with  each  component  of  the 
subset.  For  2D  digital  image  components  we  found  that  these 
structures  were  singly  connected,  so  that  by  straightforward 
application  of  the  Euler-Schlaef li  identity  we  could  compute 
the  genus  of  a  collection  of  components.  While  for  3D  digital 
components  the  structures  are  not  necessarily  singly  connected 
we  still  began  with  the  Euler-Schlaef li  equation  in  order  to 
give  insight  as  to  the  nature  of  holes  in  3D  images.  A  defi¬ 
nition  of  the  number  of  holes  in  a  subset  of  a  digital  image 
is  equivalent  to  defining  the  genus  of  the  subset  since  the 
numbers  of  objects  and  cavities  are  well  understood.  We  were 
thus  able  to  define,  and  give  an  algorithm  for  computing, 
the  genus  of  3D  digital  image  subsets. 

We  noted  in  the  introduction  that  the  2D  genus  is  not  a 
particularly  useful  property,  particularly  since  it  is  easy 
to  determine  the  number  of  objects  and  number  of  holes  inde¬ 
pendently.  In  contrast  to  the  mostly  theoretical  interest  in 
2D  genus,  the  3D  genus  is  of  more  practical  value.  Since  the 
numbers  of  objects  and  cavities  of  a  3D  image  may  be  determined 
independently,  the  genus  may  be  used  to  determine  the  number 
of  holes  in  an  object  or  objects.  Clearly  the  number  of  holes 
in  an  object  is  an  important  descriptor  that  should  be  useful 
in  many  object  recognition  tasks. 


We  have  already  noted  that  the  structures  we  have  chosen 
are  not  unique.  While  the  alternative  "structures"  mentioned 


in  Sections  3  and  5  are  not  sti LctLy  structures  as  defined 
here,  we  remind  the  reader  that  from  a  peculiar  vantage  point 
[21,  pp.  8-9]  these  too  may  be  viewed  as  structures  to  which 
the  Euler-Schlaef li  equation  is  applicable.  We  shall  not 
pursue  these  alternative  structures  further. 

Another  question  entirely  concerns  the  use  of  6-  and  26- 
connectivity  exclusively  as  complementary  types  of  connectedness. 
For  example,  we  may  define  18-connectedness  between  points 
whose  voxels  are  face-adjacent  or  edge-adjacent  [see  7] . 

Clearly  this  relation  is  an  equivalence  relation  whose  equiv¬ 
alence  classes  define  components,  and  we  can  certainly  construct 
structures  analogous  to  those  made  for  6-  and  26-components. 

We  thus  see  no  obstacle  to  the  use  of  6-  and  18-connectivity 
as  complementary  types  in  a  3D  digital  image.  We  have  chosen 
not  to  pursue  this  here  since  the  use  of  6-  and  26-connectedness 
appears  more  natural  (e.g.,  the  only  natural  distance  metrics 
on  a  3D  lattice  correspond  to  the  number  of  points  on  a  shortest 
6-  or  26-path  between  two  points  [7]).  However,  18-connectivity 
appears  to  be  equivalent  to  the  connectedness  underlying  the 
surface  topology  of  [18,19]. 

We  may  also  scrutinize  the  use  of  structures  to  model  the 
image  space.  For  example,  in  computed  tomography  a  3D  digital 
image  is  a  representation  of  a  3D  continuous  image  space  of 
objects  which  are  not  polyhedral  structures  as  used  here.  In 
Appendix  C  we  explore  the  relationship  between  G__  and  of  a 
digital  image  to  the  genus  of  a  binary  region  of  continuous 
space.  This  analysis  indicates  that  if  certain  digitization 


requirements  are  met  by  the  digital  representation  then 
and  Gg  in  fact  determine  the  genus  of  the  continuous  space,  and 
that  these  digitization  criteria  may  be  satisfied  by  using  a 
fine  enough  resolution. 

Finally,  in  Sections  4  and  5  we  introduced  a  graph  the¬ 
oretic  interpretation  of  the  genus  of  subsets  of  E,  whereby 
an  object  is  represented  by  a  graph  <•  S ,  E  (s)  >  on  the  lattice  Z 
whose  edges  are  defined  by  the  appropriate  adjacency  relation 
between  pairs  of  points  of  S.  Thus,  we  may  view  G^,  Gg , 

Gg,  G^g  as  functions  of  the  appropriate  graph. 

We  noted  earlier  that  genus  is  a  linear  set  property,  viz. 

G  (A)  +G  (B)  =G  ( Afi  B)  +G  ( AU  B) 

(see  [14]  for  a  proof  of  this).  Expanding  our  notation  to 
G  (S)  =G  ( <  S ,  E  ( S)  >  )  we  have 

G  (Aft  B)  =G  (  <  AU  B,E(AUB)  >)-G(<  A,E(A)  >  ) -G  (  <  B,  E  (B)  >  ) 
from  which  we  see  that 

G  (A(i  B)  =G  (<  0 ,  E  (AU  B)-E(A)-E(B)  >  ) 
where  0  is  the  empty  set.  This  is  similar  to  the  definitions 
in  [8]  of  simple  closed  curves  and  surfaces  by  sets  of 
ordered  pairs  (p,q)  of  points  ptS,  qfcS,  since  these  pairs 
specify  the  edges  of  E  (AU  B) -E ( A) -E (B) .  This  approach  deserves 
further  investigation. 


it**** 


Appendix  A 

The  256  binary  assignments  to  a  2*2*2  neighborhood  of 
points  in  a  3D  image  may  be  grouped  by  symmetry  into  22  dis¬ 
tinct  patterns.  The  pictorial  display  here,  taken  from  [14] 
shows  voxels  with  image  value  1  as  cubes.  In  Sections  4  and 
5  the  notation  # [n]  refers  to  the  numbering  here  of  these 
patterns.  In  patterns  15-22  only  all  hidden  voxels  have  value  1. 

(All  zeros) 

1  {=71) 


5 (=18) 


6 (=17 ) 


7 (=16) 


8 (=15 ) 


Appendix  B 

The  values  of  the  net  contribution  AG  (K.  .  )  of 

2d  x  f  3  r  k 

a  2x2x2  cube  K.  .  .  to  the  yenus  G_, (S)  of  a  subset  S  of  a 

1  i  J  »  K  2  D 

digital  image  are  tabulated  below  by  a  decimal  index  t. 

The  binary  representation  of  the  number  t  is  formed  from 
the  2x2x2  cube  of  binary  valued  points 


a - b 

g - h 


as  the  bit  string  t=abcdefgh.  The  point  with  binary  value  b 
is  considered  to  be  (i,j,k). 

Note  that  only  49  entries  for  AG26(<t>)  are  non-zero,  and 
that  -  2  *  G2g  ( <-t>)  &  1 ,  so  that  a  significant  amount  of  compression 
may  be  achieved  if  the  table  size  is  a  problem. 


AG  (<t» 


AG  (<t>) 


1 


t 

Xg  ( <t>) 

t 

AG(<t>: 

34 

0 

89 

0 

35 

-1 

90 

0 

36 

-1 

91 

0 

37 

-1 

92 

0 

38 

-1 

93 

0 

39 

-1 

94 

0 

40 

-1 

95 

0 

41 

-2 

96 

0 

42 

-1 

97 

0 

43 

-2 

98 

0 

44 

-1 

99 

0 

45 

-1 

100 

0 

46 

-1 

101 

0 

47 

-1 

102 

0 

48 

0 

103 

0 

49 

0 

104 

0 

50 

0 

105 

0 

51 

0 

106 

0 

52 

0 

107 

0 

53 

0 

108 

0 

54 

0 

109 

0 

55 

0 

110 

0 

56 

-1 

111 

0 

57 

-1 

1 1 2 

0 

58 

--1 

113 

0 

59 

-1 

114 

0 

60 

0 

115 

0 

61 

0 

1  J.  6 

0 

62 

0 

117 

0 

6  J 

0 

118 

0 

64 

0 

1 1 9 

0 

65 

0 

120 

0 

66 

0 

121 

0 

67 

0 

122 

0 

69 

0 

123 

0 

69 

0 

124 

0 

70 

0 

125 

0 

71 

0 

126 

0 

72 

0 

127 

0 

73 

0 

128 

0 

74 

0 

1.29 

-1 

75 

0 

130 

0 

76 

0 

131 

-1 

77 

0 

132 

0 

78 

0 

J  33 

0 

79 

0 

1  34 

0 

80 

0 

135 

0 

81 

0 

136 

0 

82 

0 

137 

■J 

03 

0 

1 38 

0 

84 

0 

139 

- 1 

85 

0 

.1  40 

0 

36 

0 

141 

0 

37 

0 

142 

0 

88 

0 

14  3 

0 

t 

AG  ( <t->) 

t 

laG  ( 

144 

0 

200 

0 

145 

0 

201 

0 

146 

0 

202 

0 

147 

0 

203 

0 

148 

1 

204 

0 

149 

1 

205 

0 

150 

1 

206 

0 

151 

1 

207 

0 

152 

0 

208 

0 

153 

0 

209 

0 

154 

0 

210 

0 

156 

0 

21  1 

0 

1  56 

1 

212 

0 

157 

1 

213 

0 

156 

J 

214 

0 

159 

1 

215 

0 

1  60 

0 

216 

0 

161 

-1 

217 

0 

162 

0 

218 

0 

163 

-1 

219 

0 

164 

0 

220 

0 

165 

0 

221 

0 

166 

0 

222 

0 

167 

0 

223 

0 

168 

0 

224 

0 

169 

-1 

225 

0 

1  70 

0 

226 

0 

171 

- 1 

227 

0 

1  72 

0 

228 

0 

178 

0 

229 

0 

174 

0 

230 

0 

175 

0 

231 

0 

176 

0 

232 

0 

1  77 

0 

233 

0 

178 

0 

234 

0 

179 

0 

235 

0 

180 

1 

236 

0 

181 

1 

237 

0 

182 

1 

238 

0 

188 

1 

239 

0 

184 

0 

240 

0 

185 

0 

241 

0 

186 

0 

242 

0 

187 

0 

243 

0 

188 

1 

244 

0 

189 

1 

245 

0 

190 

1 

246 

0 

191 

1 

247 

0 

192 

0 

248 

0 

193 

0 

249 

0 

194 

0 

250 

0 

195 

0 

251 

0 

196 

0 

252 

0 

197 

0 

253 

0 

198 

0 

254 

0 

199 

0 

255 

0 

Appendix  C 

In  this  appendix  we  outline  an  anai'sis  of  the  relationship 
between  the  genus  of  a  digital  image  and  the  continuous  space 
it  represents.  Let  Z*  be  a  continuous  space  composed  of 
binary  regions;  the  union  of  all  regions  with  value  1  we  call 
S* ,  and  the  union  of  the  others  S*.  We  assume  that  the  back¬ 
ground  is  a  region  of  S*. 

We  wish  to  state  the  circumstances  under  which  a  digital 
image  Z  is  a  representation  of  Z*  that  preserves  the  toplogy 
of  I*.  These  circumstances  take  the  form  of  a  local  digitiza¬ 
tion  requirement  that  we  impose  on  Z.  Let  N (p)  be  the  8- 
neighbors  of  pfcZ  including  P  when  Z  is  two-dimensional,  and 
let  N (p)  be  the  26-neighbors  of  ptZ  including  p  when  Z  is 
three-dimensional.  Associated  with  p  and  N(p)  in  Z*  we  have 
p*  and  N* (p) . 

First,  we  will  require  every  local  binary  region  of  Z* 
to  have  a  representation  in  Z : 

(i)  0(S*HN*  (p)  )  =  0(snN(p) ) 

(ii)  0(S*nN* (p)  )  =  0(SnN(p) ) 

For  the  case  where  Z  is  two-dimensional  this  is  all  we 
need.  Let  <j>  be  either  G^  or  Gg.  By  the  definition  of  <t> 
(Section  3)  it  is  not  hard  to  see  that  for  p€S 
<j>(s)=<Ms-{p})+<j>(snN(p))-<j>(  [snN(p)  ]-{p>) 


Now  every  n  point  object  S  can  be  constructed  from  an  n-1  point 
object  S— { p}  with  p  the  eastmost  of  the  northmost  points  of  S. 


We  note  that  when  S={p},  <|>(S)=1  and  G(S*)=1  by  (i)-(ii), 
and  proceed  by  induction.  By  induction  hypothesis  G((S-(p})*) 

*  <p  (S-{p) )  ,  and  since  pGS,  G  (  [S*flN  (p)  ]  -p*)  =  <f>  (  [SHN  (p)  ] -{p} )  . 

Now  G(S*HN*(p))  =  0 <S*HN*  (p)  )  -H  (S*ON*  (p)  )  and  by  (ii)  we  must 
have  H(S*nN*(p))  =  H(SnN(p))  since  H(SnN(p))  =  0.  That 
G(S*f]N*  (p)  )  =  0(S*fiN*  (p) )  =  4>  (SPIN  (p)  )  is  easily  verified  by 
using  a  program  to  enumerate  all  possible  configurations  of  an 
eastmost,  northmost  point's  neighborhood. 

The  case  where  E  is  three-dimensional  requires  an  additional 
digitization  requirement  which  ensures  that  any  hole  occurring 
in  a  small  enough  neighborhood  of  a  point  in  the  continuous 
image  exhibits  itself  in  a  well  defined  way  in  the  digital  image. 
The  form  of  this  third  requirement  is  similar  to  the  definition 
of  the  class  of  objects  with  no  holes  suggested  in  [7  ] . 

We  define  a  closed  path  TT=p.,...,p  as  one  with  p,=p  .  A 
point  p^  of  7i  is  a  path  end  if  it  has  exactly  one  neighbor 
that  is  also  a  point  of  tt,  and  is  a  path  corner  if  it  has 
exactly  two  neighbors  on  tt  which  have  a  common  neighbor  q  other 
than  tk.  Two  paths  tt=p1,  —  ,pm  and  ir'=pj,  —  ,p^  that  lie  in  C&S 
are  strongly  equivalent  in  the  component  C  if 

a)  There  exists  a  path  end  p^gir  such  that  tt  '  is  the  same 

as  tt  with  pi  and  the  repetition  of  deleted. 

b)  There  exists  a  corner  p^€n  such  that  it’  is  the  same  as 
tt  with  p^  replaced  by  q,  the  common  neighbor  of 

and  P^+2  (so  ro=n) 


'  ***’■' 


(If  occurs  several  times  on  ti  we  delete  or  divert  all  in¬ 
stances  of  it) .  We  call  the  reflexive,  symmetric,  transitive 
closure  of  strong  equivalence  equivalence  in  C.  We  say  that 
tt  is  reducible  in  C  if  it  is  equivalent  in  C  to  a  degenerate 
closed  path  consisting  of  a  single  point.  [This  definition 
of  equivalence  of  closed  paths  is  motivated  by  the  concept  of 
continuously  deformable  simple  closed  curves  in  continuous 
space] . 

We  can  now  state  our  third  requirement  as 
(iii)  If  S*HN*(p)  has  a  hole  then  there  exists  a  closed 
path  in  SflN(p)  that  is  irreducible  in  S|"|N(p) 
Intuitively,  this  means  that  if  there  is  a  simple  closed  curve 
in  S*nN*(p)  which  is  not  continuously  deformable  to  a  single 
point  without  leaving  S*HN*(p),  then  a  similar  phenomenon  should 
occur  in  SON (p) . 

The  inductive  argument  that  G(S*)  =  (j>(S)  where  <f>  is  either 
G26  or  Gg  is  similar  to  that  of  the  2D  case.  For  G2g(S)  it  is 
easy  to  see  that  since  p£S,  every  path  in  SflN(p)  is  reducible 
to  p  so  that  H(S*nN*(p)>=0,  and  that  C  (S*f|N*  (p))  =0 .  It  is  easily 
verified  (by  writing  a  program  to  test  each  configuration)  that 
for  an  eastmost,  northmost,  topmost  point  of  S,  O  (SftN  (p) )  =G_, 
(SnN(p) ) .  For  the  6-connected  case  it  is  simpler  to  show  by 
induction  that  Gg(S)-G2g(S)  =  1  (or  G2g(S)-Gg(S)  =  1). 

We  summarize  these  results  in 

Proposition  C.i.  if  s^I*  meeting  requirements  (i)  and  (ii)  if 
I  is  two-dimensional,  and  meeting  requirements  (i)-(iii)  if 


L  is  three-dimensional,  then  C.(S)=G(S*)  and  G(S)=G(S*). 


Evidently,  (i)  and  (ii)  may  be  achieved  through  the  use 
of  a  fine  enough  resolution.  While  this  is  not  as  evident 
for  (iii) ,  we  notice  that  we  could  choose  a  resolution  so 
fine  that  H (S*HN* (p))  =0  so  that  (iii)  becomes  vacuous. 

A  second  interpretation  may  be  given  to  (i)-(iii);  if  we 
define  the  genus  of  a  digital  picture  as  in  Sections  3  and  4, 
then  (i)-(iii)  define  the  classes  of  (2D  and  3D)  continuous 
image  spaces  whose  topology  is  preserved  in  their  digital  repre¬ 
sentation. 

We  leave  a  rigorous  analysis  of  this  relationship  as  an 
open  problem. 


i 
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